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Abstract 

We show that even if a lattice Dirac operator satisfies the condi- 
tions consisting of locality, free of species doublings, correct contin- 
uum behavior, 75-hermiticity and the Ginsparg- Wilson relation, it does 
not necessarily have exact zero modes in nontrivial gauge backgrounds. 
This implies that each lattice Dirac operator has its own topological 
characteristics which cannot be fixed by these conditions. The role of 
topological characteristics in the axial anomaly is derived explicitly. 
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1 Introduction 



It seems to have been commonly regarded that if a lattice Dirac operator D 
satisfies the conditions : (i) locality, (ii) free of species doublings, (iii) having 
correct continuum behavior, (iv) 75-hermiticity and (v) the Ginsparg- Wilson 
relation, then D would possess exact zero modes with definite chirality, obeying 
the Atiyah-Singer index theorem [|I| for smooth gauge backgrounds with non- 
zero integer topological charge. Explicitly, these conditions are : 

(i) D is local. 

( \D^i^{x,y)\ ~ exp(— |a; — y\/l) with / ~ a and I <^ L where L = Na is 
the size of the lattice; or D{x,y) = for \x — y\ > z with z <^ L. ) 

(ii) In the free fermion limit, D is free of species doublings. 

( The free fermion propagator D~^{p) has only one simple pole at the 
origin p = in the Brillouin zone. ) 

(iii) In the free fermion limit, D has the correct continuum behavior. 
( In the limit a ^ 0, D{p) ~ i'Jfj.Pfj, around p = 0. ) 

(iv) D is 75-hermitian. ( D"^ = 75D75. ) 

(v) D breaks the chiral symmetry according to the Ginsparg- Wilson relation 

p[. ( D75 + 75D = 2raDj5D, where r is a positive real number and a is 
the lattice spacing. ) 

However, in this paper, we show that this folklore may not be justified. 
We demonstrate that even if a lattice Dirac operator satisfies all above five 
conditions (i)-(v), it does not necessarily produce any exact zero modes in 
nontrivial gauge backgrounds. Only for topologically trivial gauge backgrounds 
with zero topological charge, these five conditions (i)-(v) can guarantee that 
the continuum axial anomaly can be recovered on the lattice. 

The basic motivation of constructing a lattice Dirac operator satisfying 
these five conditions (i)-(v) is to formulate a nonperturbatively regularized 
theory for massless Dirac fermion interacting with a background gauge field, 
as a first step towards the full theory including gauge dynamics. These five 
conditions together provide a bypass to the Nielson-Ninomiya no-go theorem 
P, and constitute the necessary conditions for D to reproduce the correct 
physics in the free fermion limit as well as in the weak coupling perturbation 
regime. However, they are not sufficient to guarantee that any D satisfying 
these five conditions must have the correct index in topologically nontrivial 
gauge fields. So far, among all D satisfying these five conditions (i)-(v), the 
Neuberger-Dirac operator is the only one which reproduces the correct index 
( as well as other important quantities such as the fermion determinant ratio 
). In the next paragraph, we show how this could be possible in principle, and 
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then present two explicit counterexamples in Section 3. Some basic properties 
of the GW Dirac operator pertinent to our discussions are collected in the 
Appendix A. 

It has been proved in ref. that for any lattice Dirac operator D satisfying 
the 75-hermiticity (iv) and the GW relation (v), the necessary condition for 
it to have nonzero index in the topologically nontrivial gauge background is 

det(I - arD) = 0. (1) 

In other words, if D does not satisfy (|I|), then the index of D must be zero. 
This condition (|1|) is equivalent to that = D{1 — raD)~^ ( which is chirally 
invariant and 75-hermitian ) has singularities. Thus, if Dc is well-defined and 
does not have any singularities, then the condition (0) is violated, and the GW 
Dirac operator D = Dc{l + arDc)^^ must have zero index, i.e., topologically 
trivial . 

If one can construct a well-defined Dc satisfying the four conditions (ii)- 
(iv) ( note that any satisfying conditions (ii) and (iii) must be nonlocal, as 
a consequence of Nielson-Ninomiya theorem ||^ ), then the GW Dirac operator 
D = Dc{l + arDc)^^ is topologically trivial and satisfies the four conditions 
(ii)-(v). Now if there exists a range of r such that D is local, then D satisfies 
all five conditions (i)-(v) but has zero index, hence, a counterexample to the 
folklore. We have at least two explicit counterexamples outlined in Section 3. 
In general, we can argue that there exists a range of r such that D is local. 
Our argument is as follows. First, we note that in the free fermion limit Dc 
has a zero mode at p = in the Brillouin zone, the condition (ii). However, 
unlike any genuine zero modes ( in nontrivial gauge fields ) which are stable 
against local fluctuations of the gauge fields, this zero mode can be elevated 
by turning on a trivial gauge background with A^{x) equal to a tiny constant. 
Henceforth we shall assume that this trivial zero mode of Dc has been elevated 
in this manner. Since both D and Dc are topologically trivial, Dc does not 
have any zero modes in nontrivial gauge fields. Thus Dc does not have any 
(non-)trivial zero modes. Then, for fixed lattice spacing a, in the limit ra ^ 1, 
D ~ (ra)~^I+ higher order corrections , a highly local operator. On the other 
hand, in the limit ra 0, D Dc, which is nonlocal. Thus, between these 
two extreme cases, there must exist a range of r such that D is local. 

Further, since the index of D = Ddl + raDc)'^ is invariant for any r [ Eq. 
( p2D ], the locality condition (i) should nothe relevant to the question whether 
D has the correct index or not. Only in the case D has the correct index ( 
i.e., the sum of the axial anomaly over all sites is correct ), then the locality of 
D comes into the play to ensure that the axial anomaly is also correct at each 
site. 

So far, we do not know of any additional constraints which can guaran- 
tee that D has the correct index satisfying the Atiyah-Singer index theorem. 
However, it is obvious that if such a criterion exists, it cannot be imposed only 
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in the perturbation regime or in the free fermion hmit hke conditions (ii) and 
(iii), since it must be pertinent to the topologically nontrivial sectors. 

Nevertheless, for any D satisfying the GW relation (v), its zero modes have 
definite chirality. So we can classify any GW Dirac operator D according to 
the ratio of its index ( n_ — n+ ) to the nonzero integer topological charge Q 
of the gauge background, where n+ ( n_ ) denotes the number of zero modes 
of positive ( negative ) chirality. The virtue of such a classification depends on 
whether this ratio is robust with respect to the topological charge Q, the size 
of the lattice ( L = Na ) and the lattice spacing ( a ) . It turns out to be the 
case provided that the number of sites in each direction is sufficiently large 
( equivalently, the lattice spacing a = L/N is sufficiently small ) such that the 
absolute value of the topological charge inside any unit cell of volume is less 
than a small number ei. 



where 



1 A. I 

P^^) = — n / , ^Vv ^7r^e^^Aatr(F^^(i/)FA^(y)) , (3) 

and F is the field tensor defined in (|^. If one uses c[D] to denote this ratio, 
then 

index(D) = n_ - n+ = c[D] Q , (4) 

where c[D] is an integer constant provided that @ is satisfied. Evidently the 
value of ei in (|^) is dependent on D. 

In topologically trivial gauge backgrounds { Q = ), D does not have any 
exact zero modes. Thus c[D] cannot be defined by Eq. (^). Nevertheless, we 
can define c[D] = 1 for the trivial sector. Then the formal definition of c[D] is 

c[D] = I ' ^ f ° (5) 

With gauge backgrounds satisfying the topological bound (0), we can clas- 
sify any GW Dirac operator according to its topological characteristics c[D] in 
nontrivial backgrounds as follows. If c[D] = 1, then D is called topologically 
proper, else if c[D] = 0, then D is called topologically trivial, otherwise D is 
called topologically improper. Evidently, we are only interested in topologically 
proper D, since this is a prerequisite for D to reproduce the continuum axial 
anomaly on the lattice. There are many examples of topologically trivial D, 
however, so far, there is only one topologically proper D which can be written 
down explicitly, namely, the Neuberger-Dirac operator Q 

D = mo{l + V), V = D^{DlDj,-^/\ < mo < 2a-\ (6) 
Dyj = —niQ + Dy/ , Dy/ : massless Wilson-Dirac operator 



a^|p(x)|<ei Vx, (2) 
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and its GW generalization 

D = 2mo(I + V)[1-V + 2mora(I + F)]"^ . 

Note that we do not have any genuine examples of topologically improper D 
except by setting mo > 2a~^ in (H) @]. 

The outline of this paper is as follows. In Section 2, we argue that if a lattice 
Dirac operator does not have the correct index on any finites lattices, then it 
cannot have the correct index in the continuum limit. Then we demonstrate 
that if D has exact zero modes in a nontrivial gauge background, then its 
topological characteristics can be revealed on a finite lattice, sometimes even 
on the smallest lattice of size 3^. We also estimate the topological bound ei 
for the Neuberger-Dirac operator. In Section 3, we present two examples of 
lattice Dirac operators which satisfy the five conditions (i)-(v), but do not 
have any exact zero modes in any nontrivial gauge backgrounds. In Section 4, 
we derive a general expression for the axial anomaly of lattice Dirac operator 
on a finite lattice, in which the role of its topological characteristics c[D] is 
displayed explicitly. In Section 5, we conclude and discuss. In Appendix A, 
we collect some basic properties of the GW Dirac operator pertinent to our 
discussions. In Appendix B, we show that nonlocal Dirac operators can have 
well-defined indices. 



2 Robustness of the exact zero modes 

It seems to have been commonly believed that a lattice Dirac operator may 
possess exact zero modes in the classical continuum limit even if it does not 
have any zero modes in nontrivial gauge backgrounds on any finite lattices. 
However, our view is that if a lattice Dirac operator ( satisfying (i)-(v) ) does 
not have the correct index on any finite lattices, then it will never have the 
correct index in the continuum limit. 

First, we set up the notations for our discussions. 

Consider a smooth gauge configuration {y4^(x)} with nonzero integer topo- 
logical charge Q in the 4-dimensional Euclidean spacetime which is temporarily 
truncated to a hypercube of size with periodic boundary conditions, i.e., 
a 4-d flat torus, T^. ( The infinite volume limit L — oo is presumably to be 
taken at the end of all calculations. ) 

The covariant differential operator acting on the matter fields is defined 

as 

T^t, = d^ + iA^ . 
The field tensor F^^,^ is defined by the commutator, 

[V^,V,] = iF^, . (7) 
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The topological charge density of the gauge background is 

P(^^ = ■:^^f^uXatT^{F^u{x)Fxa{x)) (8) 

where tr denotes the trace over the gauge group space. Then the topological 
charge 

Q = [ d'^x p{x) = [ d^x -^e^^xatT^{F^u{x)Fxa{x)) (9) 

is an integer. 

The massless Dirac operator in continuum is 



which is chirally invariant ( P75 + 751^ = ) and antihermitian ( = —V 
The Dirac matrices are defined by 




(10) 



where 



Explicitly, we choose ai,a2,o'3 to be the Pauli matrices, and (74 = il where I 
is the 2x2 unit matrix. Then the chirality operator is 

75 = 71727374 = ( - 1 ) ■ 

In a smooth background gauge field with nonzero integer topological charge 
Q, T) has zero eigenvalues and the corresponding eigenfunctions are chiral. The 
Atiyah-Singer index theorem asserts that the difference of the number of 
left-handed and right-handed zero modes is equal to the topological charge of 
the gauge configuration, 

— = Q . (13) 

The lattice regularization of the gauge background amounts to discretize 
the 4-d torus into a hypercubical lattice with sites in each direction ( i.e., 
L = Na ), and then transcribe the gauge fields {y4^(x)} to the link variables 
{U^{x)}, where U^{x) is the link variable defined on the link {x,x + afi). We 
choose the lattice spacing a sufficiently small ( i.e., is sufficiently large for 
fixed L. ) such that the variations of {A^{x)} can be tracked by the link vari- 
ables {U^{x)}. For any lattice Dirac operator D ( satisfying the five conditions 
(i)-(v) ), our primary concern is whether D can produce exact zero modes sat- 
isfying the Atiyah-Singer index theorem (p!3|) . If the index of D always agrees 
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with the topological charge Q for any gauge background satisfying the bound 
(H), then D is called topologically proper. It is evident that the value of ei 
depends on Z). A rough estimate of ei can be obtained by the following simple 
prescription. 

Consider the following gauge configuration with constant field tensor and 
topological charge (5 = 1, 



Al{x) 



Al{x) 



27rx4 



L2 







where is any one of the generators of the gauge group with normalization 
tr(t°t'^) = 5°^. The nonzero components of the field tensor are 



Fl2{x) = F34(X) 



L2 ^ - 



Then the corresponding link variables on the lattice ( = 0, 
are : 



Ui{x) 

U2ix) 

Usix) 
UAx) 



exp 
exp 
exp 
exp 



. 27rxi 



27iaxi „ 
f 



i6. 



X4„{N-l)a 



271X3 



Using the notation 

Uf,^{x) = U^{x)U^{x + afi)Ul{x + au)Ul{x) 



(14) 
(iV-l)a) 

(15) 
(16) 
(17) 
(18) 

(19) 



to stand for the path-ordered product of link variables around the plaquette 
p = (x, /t, i>), we obtain 



Ui2{x) = Uu{x) = exp 



2TTa^ 



exp 



ia'F 



(20) 



and Uis{x) = Uu{x) = U23{x) = U24,{x) = I. 

Then we test whether D has exact zero modes on finite lattices by succes- 
sively increasing the number of sites N in the sequence iV = 3,4,5,---. 

If D is topologically proper, then the exact zero mode would emerge after 
N > Ni for a certain finite integer Ni, ( i.e., for the strength of the nonzero 
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field tensor inside each plaquette, a^|Fi2| = a^|F34| = 27r/A^^, is less than 
2t[/NI ). This yields a rough estimate of ei ~ '^/{Ni + 1)"^ for the topological 
bound (|^) such that c[D] is an integer constant. A more optimal estimate of ei 
can be obtained by introducing local fluctuations which change the topological 
charge density locally but maintain the total topological charge fixed. 

On the other hand, if D is topologically trivial, then D does not have any 
exact zero modes for any finite lattice spacing a = L/N { a.^ N is increased 
successively ). Suppose that in the limit a = ( = oo ), the exact zero 
mode emerges, however, in our view, D still does not have the correct index 
in the continuum limit. Our argument is as follows. 
[ Proof ] : 

First, let the numbers of exact zero modes of D be kj^ and k_ in a nontrivial 
gauge background with topological charge Q on a finite lattice, and they do 
not yield the correct index. ( i.e., index(D) = k_ — k+ ^ Q. ) Suppose that 
the index of D remains the same for any finite lattice spacing a = L/N as N 
is increased successively, except at a = where the correct numbers of exact 
zero modes ( n± ) emerge and the index becomes n_ — n+ = Q. This implies 
that the index of D must undergo a discontinuous transition from the integer 
fc_ — fc+ to another integer n_ — n_|_ at a = 0, and correspondingly, D becomes 
Di. Thus the index(Di) = n_ — does not constitute a limiting point of the 
index(D) = k_ — k+, because a discontinuity has occurred. Since the index of 
Di is not equal to the index of D, they are in two different topological phases, 
or in other words, not in the same universality class. It follows that Di cannot 
be the continuum limit of D. This completes the proof. 

In short, if D does not have the correct index at finite lattice spacing, then 
D cannot have the correct index in the continuum limit. 

Therefore we conclude that the index of a lattice Dirac operator ( satisfying 
the conditions (i)-(v) ) is invariant for any nontrivial gauge backgrounds with 
the same topological charge provided that the topological bound (Q) is satisfied. 

Next we demonstrate that if a lattice Dirac operator is topologically proper, 
then it can produce the correct index on a finite lattice, sometimes even on a 
very small lattice. For example, on the 3^ lattice, the Neuberger-Dirac operator 
@ with mo = 1 has one exact zero mode in a gauge background with Q = 1. 
We also estimate the value of its ei in the topological bound (|^). 

For simplicity, we consider the following U{1) background gauge fields on 
the 4-d flat torus ( G [0, L], = 1, ■ ■ ■ , 4 ) : 

aM.) = -!^|^+4»)sin(5p.,) (21) 

(0) . /27mi 



aA2ix) = A^^'sini^-j-^xij (22) 
aA,ix) = -^ + 4°^sin(^.,) (23) 
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aA^ix) = Afsinf^xs) (24) 



where qi and q2 are integers. The global part of the gauge background is 
characterized by the topological charge 

Q = J d^x e^uXa F^y{x)Fxa{x) = qi q2 ■ 

The local parts are chosen to be sinusoidal fluctuations with amplitudes Ai'\ 
A^2\ Af and Af\ and frequencies ^ and ^ where rii, na, 

and n4 are integers. The discontinuity of Ai{x) ( A^^x) ) at X2 = L { 
X4 = L ) due to the global part only amounts to a gauge transformation. The 
components of the field tensor are continuous on the flat torus. To transcribe 
the background gauge field to link variables on the lattice, we take the lattice 
sites at x^ = 0, a, (A^ — l)a, where a is the lattice spacing and L = Na is 
the lattice size. Then the link variables are 

Ui{x) = exp[iaAi{x)] (25) 

iaA2[x) + ld^2,(Ar_i)a- 



U2{x) = exp 



(26) 



L 

Usix) = exp[iaA3{x)] (27) 
Uiix) = exp iaA4(x) + i(5^4,(Af-i)a ^^^^^^ (28) 

First, we consider the simplest nontrivial gauge configuration with constant 
field tensor at each plaquette, i.e., 

27rgi/iV2, {i2,u} = {l,2}; 
a'F^,{x) = -a'F,^{x) = { 2Tiq2/N\ {/i, v} = {3,4} ; (29) 

0, otherwise. 

This amounts to setting the local sinusoidal fluctuations to zero 

( Af^ = 0, 

z = l,---,4)inEqs. (^-(^. 

In Table 1, we list the number of exact zero modes of each chirality versus 
the topological charge of the simplest nontrivial gauge background, for the 
Neuberger-Dirac operator on the 3^ lattice. The Atiyah-Singer index theorem 
(|13D is satisfied exactly for Q = 1, but not for Q = 2. This yields a rough 
estimate of ei in (^ : 1/3"^ < ei < 2/3^. It is obvious that the Neuberger-Dirac 
operator is nonlocal on such a small lattice. This demonstrates that even a 
nonlocal Dirac operator can have exact zero modes satisfying the Atiyah-Singer 
index theorem. However, the axial anomaly of a nonlocal lattice Dirac operator 
must disagree with the topological charge density of the gauge background. 

Besides the zero modes, the +2 eigenmodes of the Neuberger-Dirac opera- 
tor are also chiral W\ , and the numbers of both chiralities are listed in the last 
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92 


Q = 91 92 


n+ 


?l_ 


n+ 


"2 


c[D] 


-1 


-1 


1 





1 


1 





1 


-1 


1 


-1 


1 








1 


1 


1 


1 


1 





1 


1 





1 


-2 


1 


-2 


1 








1 


1/2 


2 


-1 


-2 


1 








1 


1/2 


-2 


-1 


-2 





1 


1 





1/2 


2 


1 


2 





1 


1 





1/2 



Table 1: The exact zero modes versus the topological charge, for the 
Neuberger-Dirac operator with mg = 1, on the 3^ lattice. The Atiyah- 
Singer index theorem Q = n„ — n+ is satisfied only for \Q\ = 1. The chirality 
sum rule (pO]) is satisfied in all cases. The last column is the topological char- 
acteristics, c[D] = (n_ — nj^)/Q. The fractional value of c[D] indicates that 
topological bound is violated. 

two columns in Table 1. It has been shown that for any D satisfying condi- 
tions (iv) and (v), every zero mode of D must be accompanied by a 1/r real 
eigenmode of opposite chirality, and they satisfy the chirality sum-rule J^, 

n_ + = + n2 . (30) 

where r = 1/2 for the Neuberger-Dirac operator. In Table 1, it is evident 
that the zero modes together with the +2 chiral modes satisfy the chirality 
sum-rule (^Of) . 

Note that the exact zero modes ( as well as the +2 chiral modes ) are 
very stable under local fluctuations of the gauge background provided that the 
topological bound (0) is satisfied. ( Otherwise, they are not regarded as genuine 
zero modes of a lattice Dirac operator in a nontrivial gauge background. ) To 
demonstrate the robustness of the zero modes, as well as to obtain a better 
estimate of ei, we turn on the local sinusoidal fiuctuations in (^TJ). With 
Q = 1 on the 3^ lattice, we increase the amplitude Ai'^ gradually until the 
exact zero mode disappears. At this point, the average of a^|p(a;)| over 
all sites gives a conservative estimate of ei. In Table 2, the amplitude A^*^ 
is listed in the first column. The frequency of the sinusoidal fiuctuation is 
211 /L. The maximum and the average of a'^|p(a:)| are listed in the second and 
third columns respectively. We note that for < Ai(0) < 1.022, the index is 

— Tij^ = 1, equal to the topological charge Q = \. However, as soon as 
reaches 1.022, the exact zero mode disappears and the index becomes zero. As 
we further increase A^i^ up to 1.090, D enters into another topological phase 
with index equal to —1. Thus, at A^^ = 1.022, the average of a''|p(a;)| over 
all sites yields a conservative estimate of ei < 0.025. To be secure, we take 
ei = 0.02. For various gauge configurations on larger lattices, we have checked 
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Af 


niax(a"'|/9(a;)|) 


< a'^p{x) > 


n+ 




i4 


1^2 


c[D] 


0.000 


0.0123 


0.0123 





1 


1 





1 


0.400 


0.0185 


0.0123 





1 


1 





1 


0.600 


0.0215 


0.0164 





1 


1 





1 


0.800 


0.0246 


0.0205 





1 


1 





1 


1.000 


0.0277 


0.0245 





1 


1 





1 


1.020 


0.0279 


0.0249 





1 


1 





1 


1.022 


0.0280 


0.0250 

















1.080 


0.0289 


0.0262 

















1.085 


0.0289 


0.0263 

















1.090 


0.0290 


0.0264 


1 








1 


-1 


1.100 


0.0292 


0.0266 


1 








1 


-1 



Table 2: The exact zero modes versus the amphtude of the local sinusoidal 
fluctuations, Ai'^ in (^I]) , for the Neuberger-Dirac operator on the 3^ lattice. 
The topological charge is Q = 1. The Atiyah-Singer index theorem Q = 
— IS satisfied for < 1.022, or equivalently, a^\pix)\ < 0.028. The 
last column is the topological characteristics, c[D] = {n_ — n^)/Q. 



that the index of the Neuberger-Dirac operator with mo = 1 is always equal to 
the topological charge of the background provided that the topological bound 

a^|p(x)| < 0.02 Vx (31) 

is satisfied. 

We note that the topological bound can be satisfied if the link config- 
uration {U^{x)} fulfils 

|Re tr(I — Ui^i^{x))\ < — — ^ for all plaquettes, (32) 

where U^y{x) is the path-ordered product of link variables around the plaquette 
p = (a;,/i, i>), as defined in ([T9|) . However, the relevant parameter for the 
topological phase transition of D is the topological charge density a*^|p(a;)| 
rather than the plaquette action |Re tr(I — U^i,{x))\, since a topologically 
proper D always has c[D] = 1 for any gauge configuration satisfying (^ even 
if it violates (|32D for some plaquettes. 

Finally, we note that the topological bound (^ or (^) does not guarantee 
that D is local, but only that c[D] is an integer constant. In general, if a 
lattice Dirac operator is local in the free fermion limit, then it would maintain 
its localness in a gauge background with sufficiently small field strength at 
each plaquette, i.e., 

||2 — U^^{x)\ \ < e for all plaquettes, (33) 
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where the matrix norm ||M|| can be defined to be the square root of the 
maximum eigenvalue of M"''M, i.e., ||M|| = ^J \max{,M^ M) . Evidently the 
value of e is dependent on D. 

The locality bound (P^) implies that the topological charge inside any unit 
cell satisfies 

a^|p(a;)|<e2 Vx, (34) 
62 = ^tr(I)e^ (35) 

where tr(I) = for link variables in the fundamental representation of the 
gauge group SU{Nc). For the Neuberger-Dirac operator with mo = 1, e = 
1/6(2 + v^) ^ 0.0488 @, i, which gives ~ 5.43 x 10"^ for = 3. This 
shows that €2 is much less than ei ^ 0.02 in the topological bound. Evidently, 
for any topologically proper D, the inequality 62 < ei must hold, since the 
topological bound only requires that the index ( i.e., the sum of the axial 
anomaly over all sites ) is equal to the topological charge, while the locality 
bound further guarantees that the axial anomaly Al{x) at each site is in good 
agreement with the topological charge density p{x) of the gauge background. 

We summarize the main theme of this section as follows. Consider a fixed 
gauge background with nonzero integer topological charge Q on the 4-torus 
of size L^. Then we impose a lattice on the 4-torus with lattice spacing a, 
and the number of sites in each direction, ( i.e., L = Na ). Given any 
lattice Dirac operator D, we can test whether the index of D agrees with Q by 
successively increasing N ( equivalently, decreasing lattice spacing a = L/N ) 
in the sequence = 3, 4, ■ ■ -. If D is topologically proper, then there exists a 
finite integer A^^i such that for N>Ni{a<L/Ni), the index of D is equal to 
Q. However, when is only slightly larger than A^^i ( equivalently, a is slightly 
less than L/Ni ), D may not be local yet. In general, if D is local in the free 
fermion limit, then there exists another integer A'"2 > A^^i such that for N > N2 
{ a < L/N2 ), D is local. However, if D is topologically trivial, then D does 
not have any zero modes for any finite N, even though it has become local for 
N > N2. Whether D has any zero modes at iV = 00 ( a = ) really does not 
matter at all, since in case it has, then it cannot be the continuum limit of D. 

So far, among the lattice Dirac operators which satisfy the five conditions 
(i)-(v), the Neuberger-Dirac operator is the only one which is topologically 
proper for gauge configurations satisfying the topological bound : a^|p(x)| < 
0.02. In the next section, we present two explicit examples which satisfy these 
five conditions (i)-(v), but do not have any exact zero modes in nontrivial 
gauge backgrounds. 
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3 Topologically trivial lattice Dirac operators 



It has been proved in ref. ||^ that for any chirally symmetric and 75-hermitian 
Dirac operator Z)c, if it is well-defined ( without any singularities ), then the 
GW Dirac operator D = Dc{l+raDc)~^ has zero index in any nontrivial gauge 
background. 

Therefore, if one can construct a well-defined Dc satisfying the three con- 
ditions (ii)-(iv) ( note that any satisfying (ii)-(iii) must be nonlocal, 
as a consequence of no-go theorem ), then the GW Dirac operator D = 
Dc{l + arDc)~^ is topologically trivial and satisfies all five conditions (i)- 
(v). The locality of D in the free fermion is ensured by choosing r in the 
proper range. Then it would be local in a gauge background which satisfies 
the locality bound (^). In this section, we present two examples from recent 
studies [0, y]. 

In ref. a well-defined Dc satisfying the three conditions (ii)-(iv) has 

been constructed, 







-D 



(36) 



where 



Dl = ((T ■ ty^ w - Vw^Vl + w-^ t2 ^-1 



(37) 
(38) 



1 



ia-t){a^ -t) 



(39) 
(40) 



7/. 




and 



w[x. 



(41) 



It is evident that Dr 



is 75-hermitian since it is chirally sjTumetric and 



antihermitan. It has been shown []T0| that in the free fermion limit, Dc (|36| ) is 
free of species doublings at finite lattice spacing and agrees with 'y^dfj_ in the 
classical continuum limit. So, D^ (^) satisfies the three conditions (ii)-(iv). 

Now we examine whether D^ ( Dl ) is well-defined. Since the naive lattice 
Dirac fermion operator 7^t^ does not have any exact zero modes in nontrivial 
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gauge backgrounds, its left-handed fermion propagator (a-t)^^ is well-defined. 
Therefore, the first factor of Dl ( pTl) is well defined. Next we examine the 
second factor of We see that the second factor is well-defined except 

for the gauge configurations which give det{w) = 0. However, these exceptional 
configurations have zero measure, thus they would never be encountered in any 
practical calculations. Therefore we have confirmed that Dc (|36|) is well defined 
and satisfies the three conditions (ii)-(iv). Then the GW Dirac operator 
D = Dc{l + arDc)~^ satisfies the five conditions (i)-(v), where the locality 
of D is ensured by choosing r in the proper range. Since Dc is well defined. 



D = Dc{l + arDc 
in ref. H. Explicitly, 



is topologically trivial, according to the theorem proved 



D 



r{bb^ + r 



r 



2)-l 

2\-l 



-b{b^b + r 
r{b'^b + 'i 



2yl 

2\~1 



(42) 



where 



1 

- I w 
2 



-1 



(43) 



Another example of well-defined which satisfies the three conditions 



(ii)-(iv) is constructed in ref. 



1-1 



W 



(44) 



where is defined in (39) and W is the Wilson term |21 



W{x,y) 



x-tj.,y 



(45) 



The Dc in (|4^) is 75-hermitian, since it is chirally symmetric and antihermitian. 
It is nonlocal due to the factor (7^t^)~^ in the second term. But it is well- 
defined since the naive fermion operator 7^t^ does not have any exact zero 
modes in nontrivial gauge backgrounds. The free fermion propagator of Dc 
(^41) in momentum space is 



(46) 



sin(p^a), 



EuSin^(PM'^) and w = 2a ^ sin^(PM'^/2)- 



where 

The doubled modes are decoupled completely due to the vanishing of the 
factor t'^/{w^ + t^) at the 2'^ - 1 corners of the Brillouin zone. So the fermion 
propagator D^^{p) is free of species doublings. In the limit a — 0, Eq. ( ^61) 
gives Dc{p) = i^fiPfi around p = 0. Now we have confirmed that Dc (0) is 
well-defined, and it satisfies the conditions (ii)-(iv). 
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Then we substitute (|^) into the formula D = Dc{l + arDc) ^ to obtain a 
GW Dirac operator 



D 



where 



C 



W{a^t^)-'W 



(47) 



(4J 



Then D (^71) satisfies all five conditions (i)-(v), where the locality of D is 
ensured by choosing r in the proper range. Since Dc is antihermitian and 
well-defined, the index oi D = Dc{l + raDc) must be zero 

We have explicitly checked that the GW Dirac operators in ( ^21) and ( ^Tj) 
do not have any exact zero modes in a constant nontrivial gauge background 
with Q = 1, [ Eqs. ([T5|)-(p^ ], on different lattices with sizes ranging from 
AT = 3 to = 6. 

We note that even though these two GW Dirac operators ( ^2] ) and ( ^TD do 
not have any exact zero modes in nontrivial gauge backgrounds, the five condi- 
tions (i)-(v) do guarantee that they can reproduce the correct axial anomaly 
in trivial gauge backgrounds. This will be demonstrated explicitly in the next 
section. 



4 The axial anomaly and topological charac- 
teristics 

In this section, we derive a general expression for the axial anomaly of lat- 
tice Dirac operator D on a finite lattice, in which the role of the topological 
characteristics c[D] is displayed explicitly. 

4.1 In continuum 

First, we discuss the axial anomaly of the massless Dirac operator in a gauge 
background in the continuum, i.e., the 4-d fiat torus ( ) with size L^. The 
anomaly equation can be written as 

{d,jl{x)) = 2^(x) + 2^[0;(x)]t0^(a;) -2 5:[0i(x)]t0i(a;) (49) 

s=l t=l 

where jf^{x) = ijj{x)'y^'y5ip{x) is the axial vector current, {d^j^{x)) is the 
fermionic average of the divergence of axial vector current, 0^ and are 
normalized eigenfunctions of the zero modes with chirality +1 and —1 respec- 
tively, and A{x) is the axial anomaly. Integrating the anomaly equation ( ^9]) 
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over the 4-torus, the Lh.s. vanishes automatically, and the r.h.s. gives 

/ d'^x A(x) = n_ — . (50) 

Using the Atiyah-Singer index theorem (|13|), Eq. ( ^0]) becomes 

J^^ d^x A{x) = Q = J d^x p{x) , (51) 

where p{x) is the topological charge density defined in (H). Then (pID yields 

d^x {A{x) - p{x)) = . (52) 

According to the divergence theorem, ( |52D implies that 

A{x) = p{x) + df,g^{x) , (53) 

where the current gfj.{x) is gauge invariant and continuous on the 4-torus. 

For weak gauge backgrounds ( |74^(a;)| ^ 1/L ) with zero topological 
charge, A{x) can be evaluated using weak-coupling perturbation theory, and 
the result is equal to the topological charge density of the gauge background, 

A{x) = p{x) = -^e^^xat'c{F^u{x)Fxa{x)) . (54) 

Therefore d^g^{x) = in this case. 

Next we consider a nontrivial gauge background with constant field tensor 
Ffj,iy{x) = F° . Since field tensor is constant, A{x) must be constant. Thus 



A{x) can be solved directly from Eq. (pl|), 

^^"^^ = ^eM-A.tr(F;,F°J = po • (55) 

Again, dfj,gfj,{x) = in this case. 

It can be shown that d^g^{x) = for any smooth gauge backgrounds. 

4.2 On a finite lattice 

Next we regularize the theory by discretizing the 4-torus into a hypercubical 
lattice with sites in each direction ( i.e., L = Na ). The gauge fields {A^{x)} 
are transcribed to link variables {Un{x)} on the lattice. For any lattice Dirac 



operator D, the anomaly equation on the finite lattice is |]T2 

a'{d,4{x)) = 2a'AL{x) + 2 ^[0;(x)]t</,^(a;) - 2 ^[0*_(x)]t</,*_(x) (56) 

s=l t=l 
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where the axial anomaly is 



a^AL{x) = lim - tr \{BD-^){x,x) + {b-^B){x,x] 



(57) 



Here D = D + m, and B is an irrelevant operator denoting the chirality 
breaking of D, 



D-f5 + -i^D = B . 



(5^ 



In particular, for GW Dirac operators, B = 2raD'j^D, then (|57|) reduces to 

[0, HI, i 

a'^Aiix) = r a tr [y5D{x,x)] . (59) 

Summing the anomaly equation ( pBD over all sites of the finite lattice, then 
the l.h.s. vanishes and the r.h.s. gives 



— Ha 



X 

From (^) and (^, we obtain 

Y,a'AL{x)=c[D]Q 

X 

Now we can write 

Q = IZ"Vi(a;) , 



(60) 



(61) 



(62) 



where Pl{,x) is the topological charge density of the gauge background on the 
finite lattice. Note that we do not need an explicit expression of Pl{x), which 
is supposed to be expressed in terms of link variables. 
Then from (^) and (|62D, we obtain 

Y^iALix) - c[D]pLix)) = . (63) 

X 

The general solution of ( |63D is 

Al{x) = c[D]pl{x) + g{x) , (64) 

where g{x) satisfies J2x9{^) = 0- We can rewrite g{x) in a form analogous to 
d^g^{x) in (H), 



9{x) = J2l9,^{x) - g^{x - ap)] = d^g^{x) , 



(65) 



where g^{x) satisfies the periodic boundary conditions : g^{x + Lv) = g^{x) 
for /i, = 1, ■ ■ ■ , 4. 
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Now we write 

pLix) = p{x) + f{x) , (66) 

where p{x) is the topological charge density inside the unit cell of volume 
centered at the site x, defined in Eq. (H), and f{x) denotes the difference 
between Pl{x) and p{x). Summing Eq. (^) over all sites and using the relation 
J2x(^'^Pl{x) = J2x(^'^p{x) = Q, we obtain 

E/(^)=0- (67) 

X 

Therefore, we can rewrite f{x) in terms of a total divergence term 

fi^) = E[/m(^) - - = dMx) , (68) 

where /^(a;) satisfies the periodic boundary conditions : /^(x + LP) = //^(x) 
for /i, z/ = 1, ■ ■ ■ , 4. Substituting (|66| ) into (0), we obtain a general expression 
for Al{x), 

Al{x) = c[D]p{x) + d^G^ix) , (69) 

where 

d,G,{x) = Y.[G^{x) - G,{x - afi)] = c[D]f{x) + g{x) . 

If the variation of p{x) is very small across the unit cell centered at x ( 
which must be the case in the limit ^ oo ( a ^ ) for a smooth background 
), then p(x) ~ p{x), and Al{x) becomes 

Al{x) - c[D]^e^,xMF'^.u{x)Fx.{x)) + d^G^{x) . (70) 

where tr denotes the trace over the gauge group space. 

Equation ( pP| ) is a general expression for the axial anomaly of any lattice 
Dirac operator on a finite lattice. Note that we have not imposed any condi- 
tions on D except assuming that the exact zero modes ( if any ) of D have 
definite chirality ( this is the case for GW Dirac operators ) which has been 
used in deriving the anomaly equation (^). 

From (0), we immediately see that Al{x) can recover the continuum axial 
anomaly if and only if c[D] = 1 and d^G^{x) = for all x. However, satisfying 
the five conditions (i)-(v) does not guarantee that c[D] = 1 and dfj,G^{x) = 
in nontrivial gauge backgrounds. 



17 



Next we try to understand the total divergence term d^G^{x). Let us 
consider a nontrivial gauge configuration with constant field tensor. Then the 
topological charge density is constant and the topological charge is, 

X 

Po = ^^.uxMK-^Fl) ■ (71) 

An example is given in Eqs. ([l5|)-(p!8D. 

Since is constant, one may expect that Al{x) is constant too. How- 

ever, this is true only if D is local. In general, if D is local in the free fermion 
limit, then it would be local in a gauge background which satisfies the locality 
bound ([331) . If we assume that is the case, then Al{x) is constant and it can 
be solved directly from Eq. (§T]), 

ALix) = c[D]pLix) = c[D]po . (72) 

Comparing ( |72D to ([TOD, we obtain d^G^{x) = provided that c[D] ^ 0. If 
c\D\ = 0, then Eq. (|6l| ) becomes 

Y,Al{x) = 0, (73) 

X 

which implies that Al{x) = if Al{x) is constant. However, in a gauge 
background with nonzero constant field tensor, Al{x) cannot be zero at all 
sites, since D must at least have some responses to the gauge background 
except for D = 0. Therefore, in the case c[D] = 0, Al{x) cannot be identically 
zero. Thus ( [73D implies that Al{x) can be written as a total divergence, 

Al{x) = d,g,{x) . (74) 

Therefore, in a nontrivial gauge background with constant field tensor, the 
axial anomaly on a finite lattice is 

Thus a topologically trivial D cannot have the correct axial anomaly in non- 
trivial gauge sectors. 

Now suppose D is topologically proper. If we introduce very small local 
fiuctuations of p{x) on top of the constant po such that J2x^p{x) = 0, then 
evidently the axial anomaly of D becomes 

Al{x) = p{x) . (76) 

Further, one can deduce that Eq. (^) holds for any gauge configurations 
satisfying the locality bound (p^). 
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In the trivial sector ( Q = ), for weak gauge backgrounds ( |a74^(a;)| << 
1/N ) at finite lattice spacing or smooth backgrounds in the classical continuum 
limit, the axial anomaly Al{x) can be evaluated using perturbation theory ( 
For examples, see ref. [|l^-[^ ). If D satisfies the five conditions (i)-(v), then 
Al{x) is equal to the topological charge density of the gauge background, 

•^l{x) = -^e^uXatT{F^u{x)Fxa{x)) . (77) 

Thus even if D is topologically trivial in nontrivial gauge backgrounds, it can 
have correct axial anomaly in the trivial gauge sector. 



4.3 Examples 

In the following, we illustrate the results of ([T5|) and ( 177|) by explicit examples. 
First, for the purpose of comparison, we also consider the Wilson-Dirac lat- 



tice fermion operator ||2T| which does not satisfy the Ginsparg- Wilson relation 
(v) but the other four conditions (i)-(iv). Explicitly, the massless Wilson- 
Dirac operator can be written as 

Dw = 1,t^ + W , (78) 

where and W are defined in Eqs. (^) and (^) respectively. According to 
formula i^7\), the axial anomaly of can be written 



Awix) 



^^(4 tr[Z},^i(x,x)75] - tT[Dj{x,x + afi)-f5U^{x)] 



4 M 



- tr[D^^(x + a/t, x)75f/^(x)] - tr[D^^{x - ap,, x)'y5U^{x - afl)] 

- tT[D:^^{x,x - afi)-f5Ul{x - ap,)]). (79) 



where tr denotes the traces over the Dirac space and the gauge group space. It 
is well known that does not possess any exact zero modes in topologically 
nontrivial gauge fields, thus c[Z)vf] = 0. Therefore its axial anomaly must 
disagree with the topological charge density in nontrivial gauge backgrounds, 
according to (|70D . However, one may wonder whether can produce the cor- 
rect axial anomaly in trivial gauge backgrounds satisfying the locality bound. 
For this purpose, it suffices to consider the trivial f/(l) gauge backgrounds 
with local sinusoidal fluctuations on a 2-dimensional fiat torus of size L x L : 

aAi(x) = — ^ + 4 sin ^-^xa (80) 



L ' ' V L 
aA^ix) = — ^ + 4 ^ sin f -^si ) (81) 
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The corresponding link variables on the 2d lattice ( Xf^ = 0, a, (A^ — l)a ) 
are 



Ui{x) = exp [ iay4i(x) ] , 
U2{x) = exp [ iaA2{x) ] . 

The topological charge density on the 2d torus is 

The topological charge density inside the square centered at x is 

1 /•xi+a/2 rX2+a/2 

p{x) = — dxi / dx2 p{x) (83) 

O Jxi—a/2 Jx2—a/2 



1 



a^TT 



(0) . fnniaX f2nni \ (o) . f nn2a\ (2Txn2 
A2 sm I ^ I cos I ^ Xi j — A\ sm I ^ I cos I ^ 3:2 



M) 



In Fig. |I], we plot ^^^(2^) for each site on a 12 x 12 lattice with lattice 
spacing a = 1, comparing with the topological charge density p{x), in a trivial 
gauge background [ Eqs. ( pO| ) and ( [8T| ) ] with parameters Af"* = 0.3, ^2°'' = 0.4, 
rii = n2 = 1, /ii = 0.1 and h2 = 0.2. Note that in this case, the difference 
between p{x) and p{x) is very small, always less than one percent. The position 
of a site with coordinates (xi, X2) is represented by an integer x = I2{x2 — 1) + 
Xi, as the x-coordinate in Fig. 1. The axial anomaly Aw{x) is denoted by 
squares. The topological charge density p(x) (^) of the gauge background is 
denoted by circles. The line segments between circles are inserted only for the 
visual purpose. Evidently Aw{x) disagrees with the topological charge density 
p{x). The discrepancies are due to the presence of the fermion doublers which 
decouple completely only in the limit a — > 0. 

The deviation of the axial anomaly of a lattice Dirac operator in a gauge 
background can be measured in terms of 

1 |^.(x)-p-(x)| 
\p{x)\ 

where Ng is the total number of sites of the lattice, and p{x) is the topological 
charge density inside the unit cell of volume a"^ centered at x. 

For the axial anomaly of the Wilson-Dirac operator as shown in Fig. |l|, the 
deviation is 6w = 0.45. One expects that the deviation 6w goes to zero only 
in the limit — 00 ( or a — ). 

Next we consider the lattice Dirac operators [ (^) and (^) ] in Section 
3, which satisfy all five conditions (i)-(v), but both are topologically trivial. 
According to (^7|), the axial anomaly of (E2]) is 



Atix) = tr \{bb^ + r^)-^ - {b^b + r^)"^! (x, x) , (86) 
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and that of (P7| ) is 



Xc(x) = tr [c^Cil + r^C^Cy' - CC\l + r^CCY^] (x, x) , (87) 

where tr denotes the traces over the Dirac space and the gauge group space. 
The axial anomahes At{x) and Awdx) are plotted in Fig. ^ and Fig. |^ 
respectively, for the same gauge background as in Fig. ||. The deviation for 
the axial anomaly At{x) shown in Fig. ^ is 5 = 0.121, and for Awc{x) shown 
in Fig. 1^ is (5 = 0.068. ( For comparison, the deviation for the axial anomaly 
of the Neuberger-Dirac operator with mg = 1 is 5 = 0.049 for the same gauge 
background in Fig. |1] ). Evidently both At{x) and Awc{x) agree with p{x) ( 
p{x) ) for the trivial gauge background. They provide explicit demonstrations 
of ( |77D which holds for any D satisfying the five conditions (i)-(v), and for 
any trivial gauge backgrounds satisfying the locality bound. It is instructive to 
compare the axial anomalies of the GW Dirac operators in Fig. ^ and Fig. |^ 
to that of the Wilson-Dirac operator in Fig. |l|. The former ones can reproduce 
the continuum axial anomaly even on a finite lattice while the latter cannot. 

In Fig. ^, we plot the axial anomaly in a nontrivial gauge background ( 
Q = 2 ) with constant field strength on a 12 x 12 lattice, for the topologically 
proper Neuberger-Dirac operator, as well as the topologically improper GW 



Dirac operators (^) and (|47|) , respectively. It is clear that the axial anomaly 
of the Neuberger-Dirac operator agrees with the constant topological charge 
density po = 1/72 at each site, while those of the GW Dirac operators (^2]) 
and ( P^ ) are in complete disagreement with po- This provides an explicit 
demonstration of Eq. (|75|) for two different cases c[D] = 1 and c[D] = 0. For 
the case c[D] = 0, it shows explicitly that the axial anomaly is in the form 
d^g^{x), which oscillates around zero with its sum over all sites equal to zero. 

In passing, we note that the topological bound of the Neuberger-Dirac 
operator with mo = 1 in two dimensions is 



a^|p(x)| < ei ~ 0.28 V a; , 
while the locality bound is 

||I — < e = ~ 0.29 for all plaquettes , 

2 + v2 

which is transcribed into 

a^|p(x)| < ea ~ 0.046 M x . 
Again, €2 is much less than ei. 
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5 Summary and Discussions 



Given any lattice Dirac operator Do satisfying the four conditions (i)-(iv), 
one can construct a chirally symmetric lattice Dirac operator Dc satisfying the 
three conditions (ii)-(iv), i.e., 

= 275^0(75/^0 - Do^.r'Do , (88) 

which must be nonlocal as a consequence of the Nielson-Ninomiya theorem . 
An example of such Dc is given in Eq. (0). Then the GW Dirac operator 
D = Dc{l+raDc)~^ satisfies all five conditions (i)-(v), where the locality of D 
is ensured by choosing r in the proper range. It is obvious that one can replace 
(PSD by any D^ which satisfies the three conditions (ii)-(iv). An example is 



given in Eqs. (|36|)-(|37D. 

If Dc is well-defined ( without any singularities ) in nontrivial gauge back- 
grounds, then the index of D = Dc{l+raDc)~^ is zero 0. In other words, even 
if D satisfies all five conditions (i)-(v), D does noi necessarily have exact zero 
modes in nontrivial gauge fields. Two explicit examples have been presented in 
Section 3. At present, we do not know how to formulate additional constraints 
( in useful forms ) such that D can produce exact zero modes satisfying the 
Atiyah-Singer index theorem. Even though such conditions will be found in the 
future, it is still convenient for us to introduce the topological characteristics 
c[D] associated with each lattice Dirac operator, as defined in (^. For a given 
D, c[D] is an integer constant provided that the topological charge inside any 
unit cell is small enough, i.e., a'^|p(x)| < ei. ( For Neuberger-Dirac operator 
with mo = 1, ei ~ 0.02 ). However, in general, one does not need to go to 
the limit a = ( or = 00 ) in order to reveal the topological characteristics 
of any D. Our view is that if a lattice Dirac operator ( satisfying (i)-(v) ) 
does not have the correct index on any finite lattices, then it will never have 
the correct index in the continuum limit. An argument has been presented in 
Section 2. 

Since the index of D = Dc{l+raDc)~^ is invariant for any r, the topological 
characteristics c[D] is also invariant from r = 0, D = Dc { a. nonlocal operator 
) up to ra ^ 1, i5 ~ (ra)^^I + higher order corrections ( a highly local 
operator ). Thus the locality condition (i) does not form a constraint to c[D], 
even though it is crucial to recover the continuum axial anomaly on the lattice 
when c[D] = 1. Further, one can always transform any lattice Dirac operator 



into a GW Dirac operator with the same index [|ri[], so the GW relation (v) 
also does not form a constraint to c[D]. Therefore, we conclude that among 
these five conditions (i)-(v), only three of them [ i.e., (ii), (ill) and (iv) ] 
provide constraints to c[D]. However, satisfying all of them does not guarantee 
that c[D] = 1. Additional constraints are needed. They should be formulated 
explicitly with exact solutions rather than a set of equations with only implicit 
or approximate solutions. From this viewpoint, the Neuberger-Dirac operator 
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is much more than just a solution of the GW relation, since the former has the 
proper topological characteristics, while the latter does not necessarily have. 
This seems to indicate that the Overlap formalism |P3|, E3] ( and the Domain- 



Wall fermion |^ ) indeed plays the fundamental role in solving the problem 
of chiral fermions on the lattice. 



It is instructive to compare Liischer's definition of topological charge |22 
to that defined by the index of the Neuberger-Dirac operator. The former 
definition relies on the fact that if the gauge field satisfies the bound 

|tr(I — U^u{x))\ < e for all plaquettes, (89) 

then any link configuration ( with nonzero integer topological charge ) cannot 
deform into the trivial configuration ( with all link variables equal to the iden- 
tity ), and vice versa. Thus the gauge hnk configurations can be decomposed 
into disconnected topological sectors as in the continuum, and the topological 
charge is a well-defined integer even on a lattice. For gauge fields satisfying the 
bound (|8^) , Liischer derived an expression for the topological charge density 
on the lattice Pl(x) ( i.e., a~*q{n) of Eq. (32) in ref. |]22[ ), which satisfies 
Q = J2x(^^Pl{x), and it tends to p{x) smoothly in the continuum limit. That 
is, even on a finite lattice, Pl(x) is already in good agreement with p{x). On 
the other hand, the axial anomaly of a topologically proper D ( assuming it 
is local in the free fermion limit ) is Al{x) = p{x) ~ p{x) provided that D 
is local, i.e., the gauge background satisfies the locality bound (^). Since 
Al{x) ^ Pl{,x) for a topologically proper D [ see Eq. ( |5^ with c[D] = 1, and 
g(x) ~ ], it follows that the locality bound (|33|) of the Neuberger-Dirac opera- 
tor is compatible with Liischer's bound (^9]). However, if one is only interested 
in the topological charge, then the index of the Neuberger-Dirac operator is 
already equal to the topological charge for gauge configurations satisfying the 
topological bound (^ or ( |32D which is less restrictive than the locality bound 
(p3D . This seems to suggest that if one uses the index of the Neuberger-Dirac 
operator |Q to define the topological charge of a gauge background, then the 
lattice size can be smaller ( or the lattice spacing a can be larger ) than 
that using Liischer's definition P^ . 

A general expression (p9| ) for the axial anomaly of any lattice Dirac operator 
has been derived, in which the role of topological characteristics is displayed 
explicitly. For a topologically proper D, its axial anomaly is equal to the topo- 
logical charge density p{x) for smooth gauge background satisfying the locality 
bound. Thus, lattice QCD with the Neuberger-Dirac quarks has the correct 
axial anomaly associated with the global chiral symmetry. For chiral gauge 
theories such as the standard model, the gauge anomaly must exactly cancel 
at finite lattice spacing before one is sure the existence of a nonperturbative 
regularization of the theory. For an anomaly- free fermion multiplet, it has 



been argued that the nonabelian gauge anomaly cancellation ( in the triv- 



ial sector ) holds to all orders of an expansion in powers of the lattice spacing. 
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while the abehan gauge anomaly cancellation ^ is exact at finite lattice 
spacing. In continuum, the non-abelian gauge anomaly in 4-dimensional space 
may be obtained from the axial anomaly in 6-dimensional space ppj , pp] , pT |. 
Thus its form and normalization follow from the Atiyah-Singer index theorem 
for a certain 6-dimensional Dirac operator. Presumably, a similar analysis can 
be performed on a finite lattice with the axial anomaly in the form of ([76|) for 
a topologically proper D in a gauge background satisfying the locality bound 
(p3D. It seems that the gauge anomaly can exactly cancel at finite lattice 
spacing. We intend to return to this question in a later publication. 

To conclude, we emphasize that c[D] = 1 is a stringent requirement for any 
construction of lattice Dirac operator D in 2?7,-dimensional space. It seems 
that this requirement could not be satisfied by any approximate solutions to 
the GW relation plus other physical constraints in the 2n- dimensional space. 
Unless such 2?7,-dimensional D can possess exact zero modes and the correct 
index in nontrivial gauge backgrounds on finite lattices, otherwise its axial 
anomaly will never recover the correct result in the continuum limit. 



A Some basic properties of GW Dirac oper- 
ators 

In this appendix, we review some basic properties of the GW Dirac operator, 
which are pertinent to our discussions in this paper. Most of these properties 
can be found in refs. 0, |^ . Here we assume that D satisfies the five conditions 
(i)-(v) listed in Section 1. 

The general solution to the GW relation (v) can be written as [^] 

D = Dcil + arD^y^ = {l + arD^)-^Dc (90) 

where D^. is any chirally symmetric ( -Dc75 + 75-Dc = ) Dirac operator which 
must violate at least one of the three conditions (i)-(iii) above, according 
to the Nielson-Ninomiya no-go theorem. Now we require to satisfy (ii) 
and (iii), but violate (i) ( i.e, is nonlocal ), since (|90D can transform the 
nonlocal Dc into a local D on a finite lattice for r in the proper range, while the 
properties (ii)-(iv) are preserved. Then D satisfies all five conditions (i)-(v). 

Moreover, the zero modes and the index of D^. are invariant under the 
transformation ( PU] ) [Q. That is, a zero mode of D^. is also a zero mode of D 
and vice versa, i.e., 

D(t)± = ^ D4± = 0, 

where 75 (/>± = ±0±. Then the number of zero modes for each chirality must 
be the same for D and D^, thus the index of the local D is equal to the index 
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of the nonlocal Dc, 

n+{D,) = n+{D), n^{D,) = n^{D), (91) 

index(L),) = - n+{D^) = n_{D) - n+{D) = index(D) . (92) 

Therefore, a nonlocal Dirac operator can have well-defined index, at least for 
those obtained by the topologically invariant transformation 

Dc = D{1- arD)-^ = {l-arD)-^D , (93) 

which is the inverse transform of (|90D . From the definition of the topological 
characteristics c[D], (^), we have 

c[D] = c[D,] . 

The 75-hermiticity of D (iv) is equivalent to the 75-hermiticity of Dc, 

Dl = 75De75 ■ (94) 

Then the chiral symmetry of Dc together with its 75-hermiticity implies that 
Dc is antihermitian ( Dl = —Dc ). This last property is in agreement with the 
massless Dirac fermion operator in continuum. Then there exists one to one 
correspondence between Dc and a unitary operator V such that 

Dc = M{l + V){l-V)-\ V = {Dc- M){Dc + M)-\ (95) 

where M is a mass scale and V also satisfies the 75-hermiticity V"^ = 75^^75. 
Then the general solution (pO|) can be written as 

D = M{l + V)[{l-V) +rMa{l + V)]-^ (96) 

Due to the 75-hermiticity, the eigenvalues of D are either real or come in 
complex conjugate pairs. Furthermore, from (pO|), since r is a positive real 
number, the lower bound of real eigenvalues of D is zero, thus det{D) is real 
and nonnegative, and is amenable to Hybrid Monte Carlo simulation with any 
number of flavors of dynamical fermions. For the GW relation (v) with any r 
( not restricted to 1/2 ), the analysis in ref. [ Eqs. (24)-(41) ] goes through 
with trivial modification. The main results are : 

(a) The eigenvalues of D fall on a circle with center at l/2r, and radius l/2r, 
and have the reflection symmetry with respect to the real axis. 

{(3) The real eigenmodes ( if any ) at and 1/r have definite chirality +1 or 
-1. 

(7) The chirality of any complex eigenmodes is zero. 
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(S) Total chirality of all eigenmodes must vanish. 

Tr(75) = Es 0h50s = n+ - n_ + iV+ - iV_ = 

where n+{n-) denotes the number of zero modes of positive ( negative ) chi- 
rality, and Nj^{N^) the number of 1/r modes of positive ( negative ) chirality. 
From (5), we immediately see that any zero mode must be accompanied by a 
real 1/r mode with opposite chirality, and the index of D is 

index(D) =n_-n+ = -{N_- N+) . (97) 

Now the central problem is to construct the chirally symmetric which is 
nonlocal, and satisfies (iii), (iv), and (^41). Furthermore we also require that 
Dc is topologically proper ( i.e., satisfying the Atiyah-Singer index theorem ) 
for any smooth gauge background satisfying the topological bound (0). These 
constitute the necessary requirements [jl6l for Dc to enter (|90D such that D 
could provide a nonperturbative regularization for a massless Dirac fermion 
interacting with a background gauge field. Explicitly, these necessary require- 
ments are : 

(a) Dc is antihermitian ( hence 75-hermitian ) and it agrees with '^^{dfj_+igA^j) 

in the classical continuum limit. 

(b) Dc is free of species doubling. 

(c) Dc is nonlocal. 

(d) Dc is well defined in topologically trivial background gauge field. 

(e) Dc has zero modes as well as simple poles in topologically non-trivial 

background gauge fields ( each zero mode of Dc must be accompanied 
by a simple pole of Z^c )• Furthermore, the zero modes of Dc satisfy the 
Atiyah-Singer index theorem for any smooth gauge background satisfying 
the bound (0). 

The general solution of Dc satisfying these requirements had been inves- 
tigated in ref. [[l^]. However, in general, given any lattice Dirac operator D, 



there exists a transformation T{Rc) for D such that the transformed Dirac 



operator Dc = T{Rc)[D] is chirally symmetric [12, 11 1. 



B Indices of nonlocal Dirac operators 

In this appendix, we note that the locality condition (i) may not be relevant 
to the exact zero modes of a Dirac operator. Consider a massless fermion in a 
background gauge field, the Dirac operator in continuum is 

T> = l^{d^ + lA^) , (98) 
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which is local, chirally invariant ( 1^75 + 75P = ) and antihermitian ( = 
—T) ). Now we can define a new Dirac operator 

V = V{1 - rhV)-^ = {1 - rhVy^V , (99) 

where r is a positive real number and h denotes a regulator with mass dimension 
M~^. One can always choose the value of r such that V is nonlocal. Now it 
is evident that any zero mode of T) must be a zero mode of P', and vice versa, 
i.e., 

V(l)± = ^ = 0, 

where 75 0± = ±0±. Then the number of zero modes for each chirality must be 
the same for T) and "D', thus the index of the local Dirac operator T) is equal 
to the index of the nonlocal Dirac operator "D', 

n+{V') = n+{V), n.{V') = n_{V) 

index(P') = n_{V') - n+{V') = n_{V) - n+{V) = index(P) . (100) 

Therefore, the nonlocal Dirac operator V ( p9D has well-defined index, and it 
also satisfies the Ginsparg- Wilson relation 

V'-f5 + = -2rbV'^^V' . (101) 

It is interesting to note that the chiral limit ( i.e., r = ) of the nonlocal V 
(p9D is the local operator T> in (0), while on the lattice, the chiral limit of a 
local GW Dirac operator D satisfying (v) is 

D^ = D{l-arD)-^ = {l-arD)-^D , (102) 

which must be nonlocal if is free of species doubling and has the correct 
continuum behavior, according to the Nielson-Ninomiya no-go theorem. Again, 
the index of Dc is equal to the index of D. Therefore, no matter in continuum 
or on a lattice, a nonlocal Dirac operator can have well-defined index, at least 
for those obtained by the topologically invariant transformations, (^) and 
( |102[ ). In other words, if a nonlocal lattice Dirac operator does not have exact 
zero modes in topologically nontrivial gauge backgrounds, then the cause may 
not be due to its nonlocalness, but its topological characteristics. 
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p(x) : topological charge density 
Ajx) : axial anomaly of Wilson-Dirac fermion 
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Figure 1: The axial anomaly Aw{x) [ Eq. d79| ) ] of the massless Wilson-Dirac 
operator [ Eq. ([78|) ] in a trivial gauge background on a 12 x 12 lattice. The 
background U{1) gauge field [ Eqs. (|80|)-(|8TD ] is specified by the parameters 



hi = 0.1, /i2 = 0.2, A^^ = 0.3, A^2^ = 0.4 and rii = n2 = 1. The axial anomaly 
of the Wilson-Dirac operator, Aw{x), is denoted by squares. The topological 



charge density p{x) [ Eq. (|83|) ] of the gauge background is denoted by circles. 



The line segments between circles are inserted only for the visual purpose. 
Evidently Aw{x) disagrees with p{x). 
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p(x) : topological charge density 

A,(x) : axial anomaly of a GW Dirac operator 
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Figure 2: The axial anomaly At{x) [ Eq. ( p6D ] of the massless GW Dirac 
operator D [ Eq. (||) ] in a trivial gauge background on a 12 x 12 lattice. The 
value of r in Z) has been set to 0.4, and there is no significant changes to At{x) 
for any r in the range ~ 0.2 to ~ 0.8. The background U{1) gauge field is 
the same as that in Fig. 1. The axial anomaly At{x) is denoted by triangles. 
The topological charge density p{x) [ Eq. ( pS] ) ] of the gauge background is 
denoted by circles. The line segments between circles are inserted only for the 
visual purpose. 
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p(x) : topological charge density 
V A^^(,(x) : axial anomaly of a GW Dirac operator 
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Figure 3: The axial anomaly Awc{x) [ Eq. [ p7\ ) ] of the massless GW Dirac 
operator D [ Eq. (|47|) ] in a trivial gauge background on a 12 x 12 lattice. The 
value of r in D has been set to 0.5, and there is no significant changes to Awc{x) 
for any r in the range ~ 0.2 to ~ 0.8. The background U{1) gauge field is 
the same as that in Fig. 1. The axial anomaly Awc{x) is denoted by triangles. 
The topological charge density p{x) [ Eq. ( ^3] ) ] of the gauge background is 
denoted by circles. The line segments between circles are inserted only for the 
visual purpose. 
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Figure 4: The axial anomalies of GW Dirac operators in a nontrivial gauge 
background with constant field tensor F12 = 7r/36 ( topological charge Q = 2 
) on a 12 X 12 lattice. The topological charge density ( p{x) = 1/72 ) of 
the gauge background is denoted by the horizontal line. The axial anomaly 
An{,x) of the Neuberger- Dirac operator agrees with 1/72 excellently, lying on 
the horizontal line. The axial anomaly A.wc{x) of ( ^Tj ) is denoted by diamonds, 
while At{x) of p^) is denoted by triangles. Evidently Awc{x) and At{x) both 
disagree with the topological charge density p(x) = 1/72. 
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